epl draft 



Can dephasing generate non-local spin correla- 
tions? 



Colin Benjamin 1 

1 National Institute of Science Education and Research, IOP Campus, Sachivalaya Marg, 
Bhubaneswar 751005, India 



PACS 74.50.+r- 
PACS 74 . 78 . Na - 
PACS 85.25.-j - 

Abstract -By examining the full counting statistics of a non adiabatic pure spin pump with 
particular emphasis on the second and third moments, it is shown that incoherent or sequential 
transport, in contrast to coherent transport, can change non-local spin shot noise cross-correlations 
from being anti-correlated to being completely correlated, a truly counterintuitive result. The third 
moment on the other hand is shown to be much more resilient and its nature remains unaltered 
in incoherent transport regime. However, phenomenologically including dephasing modifies this 
picture as both Shot noise and more so the third moment are non-trivially affected. In fact non- 
local spin correlations are completely positive for maximal dephasing. 



Introduction. — Non-local shot noise correlations, the second moment, in solid state 
devices have been studied for a long time. Some of these studies include normal metal- 
superconductor hybrid structures [I], coulomb blockaded quantum dots [2], exploiting the 
Rashba scattering [3], etc. However, an experimental demonstration has thus far been lack- 
ing. This is mainly due to the difficulty in controlling environmental effects like dephasing 
or decoherence. The origin of these environmental effects can be traced to magnetic im- 
purities in the experimental system which affects electronic spin and to temperature which 
leads to electron-phonon interaction and dephasing. It begs the question how to deal with 
decoherence and reduce it. In this work a novel scheme is proposed in which the dephas- 
ing present in such systems can be used as a resource. We particularly concentrate on 
the electronic spin. The reason for dwelling on the spin instead of charge is because there 
have been many works on the charge counting statistics however works on the full counting 
statistics for spin are less visible. However, they have been attempted in different context 
to that which is the topic of this letter. For example, in [4], the FCS of spin currents was 
first attempted, the FCS of spin transfer through ultra small quantum dots in context of 
Kondo effect was attempted in [5] while in [6] a study of FCS in interacting quantum dots 
attached to ferromagnetic leads revealed super-poissonian transport. Many works revolve 
around the non-local spin shot-noise correlations. Among the notable works on non-local 
spin shot-noise correlations mention may be made of: spin current shot noise of (i) a single 
quantum dot coupled to an optical micro-cavity and a quantized cavity field [7] , (ii) a realis- 
tic superconductor-quantum dot entanglcr [8] , and (iii) a spin transistor f9 |llOJ . In this letter 
the properties of the third moment are also calculated. The reason for looking at the third 
moment is because, for charge in contrast to spin transport, the third moment is predicted 
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to be much more resilient to decoherence . In our work we see that this statement does 
not strictly hold for third moment non-local spin correlations especially when decoherence 
is large. 

Charge or spin transport is a statistical process involving electrons carrying definite 
amounts of spin or charge, since charge or spin current fluctuates in time. Therefore, in 
addition to knowing the mean charge or spin current passing through a normal conductor 
one needs to know the noise as well as the other transport moments in order to fully char- 
acterize charge or spin electron motion. To do this one takes recourse to the full counting 
statistics(FCS), which gives us the complete knowledge about all the moments of the dis- 
tribution of the number of transferred charges or spins. The full counting statistics for a 
non-adiabatic pure spin pump is analyzed in both the completely coherent and incoherent 
transport regimes. 

In this letter we find that in the coherent transport regime the current and non-local 
spin shot- noise correlations arc similar to that in Rcf. [12] . In the sequential or incoherent 
transport regimes the odd moments are almost unchanged. In contrast the second moment, 
i.e., spin shot noise becomes completely positive. An extremely counter-intuitive result. 
In absence of any dephasing the third moment spin auto or non-local correlations do not 
change much from the coherent and incoherent transport regimes. To connect the coherent 
and incoherent transport regimes we introduce a phenomenological model of dephasing. 
In this work it is shown that for maximal rate of dephasing the coherent and sequential 
transport regimes match exactly. 

The main body of this letter starts with an explanation of the model. The coherent den- 
sity matrix equation is then analyzed separate from the incoherent density matrix equation 
to bring out the differences. Lastly we bring out a perspective on future endeavors. 




Fig. 1: The model system. A single quantum dot is connected to two leads. A large coulomb 
repulsion is assumed for the quantum dot leading to prevention of double occupation. A large 
oscillating magnetic field with strength denoted by R r f — g±HBB r f /2 establishes pure spin currents 
in the leads by pumping spins from lower to higher energy level. 

Model. — The model of Ref. [H] is the starting point. It is depicted in Fig.l. It is a 
single quantum dot connected to two leads. The single electron levels in the dot are split 
by an external magnetic field B. Thus, e-j- — = g z fisB = A(Zccman energy), where g z is 
effective electron g-factor in z-direction and [Ib is Bohr magneton. No bias voltage is applied 
across the leads. An additional oscillating magnetic field B r f(t) = (B r f cos(wi), B r f sin(wt)) 
applied perpendicularly to constant field B with frequency uj nearly equal to A can pump 
the electron to higher level where its spin is flipped, then the spin down electron can tunnel 
out of the leads. Coulomb interaction in the quantum dot is considered to be strong enough 
to prohibit double occupation. No extra electrons can enter the quantum dot before the 
spin-down electron exits. The Hamiltonian of ESR induced spin battery under consideration 
is written as: 
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+ Y.^ c U^ + h.c) + H rf (t) (1) 

77, ft, (7 

In the above equation, <£ k a (c v ,k.a) and c da (cda-) are the creation and annihilation operators 
for electrons with momentum k, spin a and energy t lh k,<? hi lead rj(= L, R) and for spin a 
electron on the quantum dot. The third term describes coulomb interaction among electrons 
on the quantum dot. The fourth term describes tunnel coupling between quantum dot and 
reservoirs. H r f(t) describes the coupling between the spin states due to the rotating field 
B r f(t) and can be written in rotating wave approximation as: 

H rf (t) = R rf (c^c di e^ + create-***) (2) 

with, ESR rabi frequency R r f = g±fJ.BB r f /2, with g-factor g± and amplitude of rf field B r f. 

The quantum rate equations for the density matrix can be easily derived as in Ref. [12] . 
poo and p arj describe occupation probability in QD being respectively unoccupied and spin- 
<j states and off-diagonal term p^\.(\.-\) denotes coherent superposition of two coupled spin 
states in quantum dot. The doubly occupied is prohibited due to infinite coulomb interaction 
U — > oo. To derive the density matrix, we proceed as follows. The time dependence can be 
removed from Eqs. [TJH], by using the following unitary transformation [13] : 

JJ — e ~^[^ ktr ,( C l,k,l C vA,±~ C l,k^ C r,,k,t) + ( C d,l C d,±- C d l1: c d,t')] ^ 

The Hamiltonian is then redefined in the rotating reference form as follows: 

, dU- 1 
H RF = V~ X HV + i——U 

dt 

+ ( v n c ! h k,* c d<? + h.c.) + R r f(c dt c di + 4 ; c dt ) 

rf,k,(T 

(4) 

In the above equation, e-f = eo — y + T' ano - = e D + — -f' wnue ^fet = (r i k + T anc ^ 

- w_ 

To get the density matrix from the above Hamiltonian, the following procedure is used. 
An electron operator affecting only the electron on the dot can be written in terms of 
\p >< p\,p = 0,til- Writing, for the annihilation operator of the dot c da = |0 >< a\, 
and for the creation operator for the dot c\ la = \a >< 0|, the Hamiltonian is rewritten in 
terms of the three states: |0 >, | f>, | 4_>, corresponding to empty state, a single electron 
with spin-up and single electron with spin-down. The doubly occupied state in the dot is 
prohibited by the fact that U is taken to be extremely large. This is the infinite-U Anderson 
model. We invoke this approximation to avoid double occupancy of the dot. This entails 
that the probability for an electron entering the dot depends on whether dot is already 
occupied or not and not on the mean occupation of the dot. Thus the Hamiltonian reduces 
to: 

H = 

T},k,<7 CT 

+ E ( V A,k,a\° >< ff l + h - C -) + R rA\ t><! I 

+ U><tl) (5) 

The elements of the density matrix p mn in dot spin basis are expectation values of operators 
\n >< m\, with n,m = 0,t,i, so we can write-poo =< |0 >< 0| >,p aa =< \a >< a\ > 
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, p aS =< \a >< a >. The time evolution of the density matrix elements can be expressed 
in terms of expectation values for new operators [17) . For instance, 

ipoo = ij t < |0 >< 0| >=< [|0 >< 0\,H] > 



Poo = i[H\0 >< 0| - |0 >< 0\H] 



i[V v *\a><0\c vk<T ~V v \0><a\cl 



= {V*G< ka (t,t)-V v G<^ vka (t,t)} (6) 

To derive the Greens functions for the dot-lead system we assume that the dot-lead 
system is weakly coupled. This means the dot Greens function in presence of tunneling 
can be written is the same form as the decoupled dot Greens function. This is also called 
the Markov approximation in which the probability of a tunneling event at a given time 
depends only on the occupation at that particular time, i.e. there is no memory structure in 
the system. Invoking Markov approximation is reasonable as tunneling happens rarely and 
the system is in same state at each tunneling event, i.e., when system is weakly coupled. 
The approximated current Green's functions are (using Ref. |17j ) as a guide we have: 

G< ka > ,0a (M') = J dh^itMW^G^^hX) 

+ g^A^m^Gl.^t 1 )] (7) 

The Goo-cr''s are the green functions for the dot, while g^ka is the Green's function for the 
77-lead in absence of tunneling. 

From the convolution theorem for Fourier transforms, 

dhA(t - - t) = J duA{u)B{-u) = J ^-A(w)B(w). (8) 

Inserting the approximated current Green's functions from EqsJT] into Eq|5] and Fourier 
transforming one gets: 

Poo - \V,\ 2 [G< a a {w){g* ka {w) - g^H) + S^M^-M - G^(w))} (9) 

The general property for Green's functions G > — G < = G H — G A , is then used- 

poo = \v v \ 2 [G^Aw)(g>M - + < fe .M(Go^M - G&„(tu))] (10) 

The lesser Green's function then becomes- 

<fc CT (*) =< c lka c nka(t) >= ie- 1 ^ 1 < cl ika c vka >= ic- <e '*'*/ 1? (e, fc<7 ), (11) 
where, /(e) is the Fermi function. Performing a fourier transformation yields 
9^ka( w ) = "2Kif v (e vk <j)5{w - e I)fe(J ),and similarly g> ka (w) = -2iri[l - f v (e vk<T )]S(w - e, lkr7 ) (12) 

Substituting the above expressions in EqsJBl and using the coupling parameter rjj(e) = 
27r Efc \V v \ 2 S(e - e vka ) gives- 
Poo = g J dw^imi - f n H)G< aa M + m,{w)G> aa {w)} (13) 
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The lesser and greater Greens functions for the dot can be derived using the same formalism 
as in Ref. [T7]. Thus, GQ aa (w) = 2mp a(J S{w — e a ), and G^ crcr (u>) = —2itip m 5(w — e a ). After 
substituting these expressions in EqfT51 and integrating gives- 

poo = £r2[(l - f v (e*))p„ a - fr,(ea)poo] (14) 

at] 

Now in Ref. [T2] the Fermi functions for the left and right leads with respect to the electron 
spin /i(e-|-) = /i?(ef) = 1 and /l(e^) = = 0. Further the coupling parameters T's 

are independent of energy which implies that the density of states and tunneling matrix 
elements are constant. This approximation is called the wide band limit. It is also assumed 
that occupations are constant in time as we are only interested in steady state result where 
this approximation is valid. The weak coupling assumption as invoked above also implies 
that no broadening of energy level occurs in the dot due to tunneling and this means that 
coupling is much smaller than temperature. Thus, 

poo = -(Tf + r«)p o + (Tf + Tf) Pu (15) 

Proceeding in exactly the same way, and using the Ref. |17j as a guide one can derive the 
other rate equations as written below. To model incoherence we turn to Ref. [15] and use 
that as a model. 

Results. — We introduce density matrices p a b(t) meaning quantum dot is on the elec- 
tronic state \a > (a = b = 0,t, I) or on a quantum superposition state (a ^ b) at time t. 
We introduce counting fields [33], Xn,o, V = L/R and a =f / \. to describe transitions from 
the dot to leads. 

Coherent transport regime:. We first deal with the coherent regime: p{t) = (poo, Pff, Pa, ^(Pti)} ^(Pt-l)) = l 
Mp{t), with 

/ -(r it + r flT ) o (r^e^ + r^e^) o o \ 

(r Lt e- l XL t +T m e~ 1 ^) -2R rf 

M= -(T Ll + T Ri ) 2R rf 

-(TLi + r Ri ) -Sesr 

\ R r f —Rrf Sesr —(^li + ^ri) } 

(16) 

and Sesr — A — uj. The normalization relation poo + l~2 aa Poo = 1 holds for the conservation 
and r., ;cr = 2ir^2 k |V^| 2 5(w — e-qko)- We assume the spin relaxation time of an excited spin 
state into the thermal equilibrium to be very large. 

We calculate the eigenvalues of Eq|T6] The minimal of these eigenvalues defines the full 
counting statistics (as, Xn<* — >■ 0, 77 = L,R;a =t)4-)- After finding this eigenvalue Evq, and 
then by using the approach pioneered in Ref. [TJ], We calculate the first, second and higher 
cumulants. Note that the approach of Ref. [T3] has been generalized in Refs. jTSJUB] to 
include both coherent and incoherent transport regimes. 

The first cumulant is defined as the current, we calculate the individual spin polarized 
currents as follows: I va = | Xr)g _i.o- The spin current is thus = — I v i, while the 
charge current is K = I^ + I^i- The second cumulant defines the shot-noise. The shot noise 
local and non-local correlations can be calculated as follows. The spin shot noise local and 
non-local correlation is what we concentrate on. S S LL = S L \ + S L ^ L — S L ^ L — S LL and S LR = 

Slr + Slr-Slr-Slr herein, S$ = 9x fj£ l<r , I*,, W -*>- Similarly the third moment 

spin correlations are calculated as follows: C* ,„„ = C^tL, + C*iL, + C£bL + C, U1 



\ r)r)' r\" 



rjrj 7] TjTj' rj" ' rjrj'rj" ' rm'r]" 1 7]7]'r}" 

/-ftJ-t i rtlYt I \ ,, r v 10T . ; ri (~<oo'o" _ d 3 Ev 

V-i" + Si'i" + ^vv'v")> wnerem ^ m < n » - dx^d XWa >dxw'a" \xv*,x n '*'<Xv"-"- 



The second and third spin cumulants are sum of the individual correlations and arc given 
as Cf = Si L + S RR + 2Sl R and CJ = C S LLL + C RRR + 3C£ Li? + 3C£ Rfl . The spin Fano 
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Fig. 2: (Color online) A comparison of coherent and incoherent transport regimes. The odd 
moments- spin currents(left) and third moment cross and auto-correlations are plotted as func- 
tion of strength of rotating field for various dephasing rates F^/r = O(black), 10(red), lOO(blue). 
Solid lines are for Coherent regime and dashed lines are for sequential transport regime. The 
parameters are F = 1, Sesr = 0. 



factor and normalized skewness also can provide more information as to how spin transport 
is affected by dephasing these are defined as- spin Fano factor C^jV and normalized spin 
skewness as Cf/P. 

Sequential or Incoherent transport regime:. To go into the incoherent or sequential 
transport regime as exemplified in Refs. [15], we use the complete coherent matrix, EqlTBl 
The coefficient matrix for incoherent transport can be obtained from Eqll6[ via setting 
$t(pfl) = and S(/9|j.) = and then solving the two simultaneous equations for $t(pt±) and 
S(p^) as in Refs. [FolfTB] . This leads to a 3X3 matrix: p(t) = (poOi Att> Pu) = ^P(*) 

/ -(r Lt + r fit ) o r^ + r^e^ \ 

M = r £t e-^t + Tate-**** -z z , (17) 

V z -z-(T Li + T Ri ) J 

and, z = ^^^r^^r 1 ^ >2 ■ The minimal eigenvalue of this equation is again what we require. 

Model for Decoherence. In order to understand the coherent and sequential transport 
regimes better and how could the transition between these two regimes be connected we 
introduce a phenomenological model of decoherence via a charge detector. In this model 
the off-diagonal elements of Liouville equation are considered to exponentially decay to zero 
with the rate or 1/T 2 , where T 2 is the spin decoherence time, i.e., in the last two rows of 
the coefficient matrix (|TB)) . T is replaced by T + r^. This method of introducing decoherence 
can be substantiated by the insertion of a quantum point contact close to the quantum dot. 
Whenever an electron enters the quantum dot the transmission through it changes. This 
charge detection leads to exponential damping of the off-diagonals as derived in Refs. [T81H9] . 
Similarly, in the sequential transport regime the decoherence factor can be introduced via 
the replacement r — > T + in z of the coefficient matrix for sequential tunneling, Eq. (|17[). 

In the incoherent regime too the spin current, spin shot- noise auto and cross correlations 
are calculated and finally the third moment auto and cross-correlations. In Fig. 2, the 
odd moments are plotted- pure spin current I s and the third moment auto C S LLL , and cross- 
correlations G S LLR . In Fig. 3, the second moment, shot noise auto 5£ L and cross-correlations 
S S LR . Finally, in Fig. 4 the spin cumulants (both second and third) and the Fano factor and 
normalized spin skewness are plotted. In the insets of the right hand panels of Fig. 4 the 
large R r j is shown. One thing which is quite clear in all figures is that dephasing washes out 
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Fig. 3: The second moment: The non-local spin shot-noise correlations plotted as function of 
strength of rotating field for various dephasing rates r^/T = O(black), 10(red), lOO(blue). Solid 
lines are for Coherent regime and dashed lines are for sequential transport regime. In the inset spin 
auto (or, local) correaltions are depicted. Parameters are F = 1,Sesr = 0. 

all features in cumulants making them just increase monotonically with increasing strength 
of the rotating field. Both the spin cumulants are zero for small R r f with increasing R r f 
the second and third spin cumulants differ markedly, with the third showing much more 
feature including a pronounced dip at just above R r f = 0.5T. The spin Fano factor plotted 
on top right hand corner of Fig 4 saturates in the large R r f limit at just above 2 spin flux 
quantums this is in agreement with the normalized skewness which also saturates just under 
this value. 

In all of these figures the results for the coherent and incoherent transport regimes are 
contrasted as dephasing rate is increased from T^/T = to 100. In both regimes the charge 
current is absolutely zero. Thus there is a pure spin current. The physics behind the pure 
spin current can be outlined as follows. In the model (Fig. 1) coulomb interaction in the 
quantum dot is strong enough to prohibit the double occupation, no more electrons can 
enter the quantum dot before the spin-down electron exits. As a result, the number of 
electrons exiting from the quantum dot is equal to that of electrons entering the quantum 
dot; namely, the charge currents exactly cancel out each other implying zero charge current. 

Conclusions. — The pure spin current obtained in our set-up is shown to give rise to 
completely positive shot noise cross-correlations in the sequential transport regime. This 
is seen to be sustained when dephasing is included. In fact maximum dephasing gives rise 
to completely correlated non-local spin correlations. An analysis of the case when there is 
no dephasing is presented in Table 1. For maximal dephasing the coherent and sequential 
transport regimes merge. 

The main result of our work is depicted in Fig. 3, this is perhaps the first work where it is 
shown explicitly that the shot noise cross-correlations turn completely positive for maximal 
dephasing (blue lines in Fig. 3) in either regime. What are the reasons for the completely 
positive shot noise cross-correlations? One can see from the formula for the spin shot noise 
cross-correlations it is a difference between same spin and opposite spin correlations. For 
maximal dephasing one notices that the magnitude of the same spin correlations, which are 
negative, is always less than that of the opposite spin case. In Fig. 2, the third moment 
auto and cross-correlations are also plotted. The odd moment doesn't change markedly as 
compared to the second moment in case of sequential transport. Dephasing however has a 
dramatic impact. The third moment auto-correlations arc completely negative as expected 
since the possibility of detecting three electrons is prohibited via Paulli exclusion while 
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Fig. 4: Spin cumulants: C|(top left panel) and Cf (bottom left panel) alongwith their respective 
Normalized counterparts, Spin Fano factor (top right panel) and Normalized skewness (bottom 
right panel) plotted as function of strength of rotating field for various dephasing rates T<p/F — 
O(black), 10(red), lOO(blue). Solid lines are for Coherent regime and dashed lines are for sequential 
transport regime, in the insets the large field limit is shown. Parameters are T = 1, Sesr = 0. 



third moment cross-correlations turn negative for maximal dephasing. We have compared 
and contrasted the absolutely incoherent and absolutely coherent regimes. An effective 
parameter which shows the transition from completely coherent to completely incoherent 
can be introduced in the coherent density matrix, Eq. 16, to model this. Phenomenologically 
introducing a spin relaxation time into the coherent density matrix does indeed show the 
transition between completely coherent and incoherent regimes attesting our results. In 
this article for the first time the dramatic nature of Non-local spin shot noise correlations 
as function of dephasing is shown. Future endeavors on effects of incoherence on different 
geometries especially including superconductors are contemplated. 



Table 1: Comparing first three moments in coherent and incoherent (or sequential )transport regimes 
for zero dephasing 



Moment 


Coherent 


Incoherent 




Pure spin current 


Pure spin current 




Non-local shot-noise 


Non-local shot-noise 




correlations positive for certain 


correlations 




range of parameters 


always positive. 


3 rd 


Third moment finite 


Third moment finite 
No qualitative change 
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